Lagrange inversion theorem

In mathematical analysis, the Lagrange inversion theorem, also known as the Lagrange–Bürmann formula, gives the Taylor series expansion of the inverse function of an analytic function.

Contents

Theorem statement

Suppose the dependence between the variables w and z is implicitly defined by an equation of the form

f(w) = z\,

where f is analytic at a point a and f '(a) ≠ 0. Then it is possible to invert or solve the equation for w:

w = g(z)\,

where g is analytic at the point b = f(a). This is also called reversion of series.

The series expansion of g is given by[1]

g(z) = a %2B \sum_{n=1}^{\infty} \left( \lim_{w \to a}\left( \frac{\mathrm{d}^{\,n-1}}{\mathrm{d}w^{\,n-1}} \left( \frac{w-a}{f(w) - b} \right)^n\right) {\frac{(z - b)^n}{n!}} \right).

The formula is also valid for formal power series and can be generalized in various ways. It can be formulated for functions of several variables, it can be extended to provide a ready formula for F(g(z)) for any analytic function F, and it can be generalized to the case f '(a) = 0, where the inverse g is a multivalued function.

The theorem was proved by Lagrange[2] and generalized by Hans Heinrich Bürmann,[3][4][5] both in the late 18th century. There is a straightforward derivation using complex analysis and contour integration; the complex formal power series version is clearly a consequence of knowing the formula for polynomials, so the theory of analytic functions may be applied. Actually, the machinery from analytic function theory enters only in a formal way in this proof, in that what is really needed is just some property of the formal residue, and a more direct formal proof is available.

Applications

Lagrange–Bürmann formula

There is a special case of Lagrange inversion theorem that is used in combinatorics and applies when f(w)=w/\phi(w) and \phi(0)\ne 0. Take a=0 to obtain b=f(0)=0. We have

g(z) = \sum_{n=1}^{\infty} \left( \lim_{w \to 0} \left( \frac {\mathrm{d}^{n-1}}{\mathrm{d}w^{n-1}} \left( \frac{w}{w/\phi(w)} \right)^n \right) \frac{z^n}{n!} \right)
= \sum_{n=1}^{\infty} \frac{1}{n} \left( \frac{1}{(n-1)!} \lim_{w \to 0} \left( \frac{\mathrm{d}^{n-1}}{\mathrm{d}w^{n-1}} \phi(w)^n \right) \right) z^n,

which can be written alternatively as

[z^n] g(z) = \frac{1}{n} [w^{n-1}] \phi(w)^n,

where [w^r] is an operator which extracts the coefficient of w^r in what follows it.

A useful generalization of the formula is known as the Lagrange–Bürmann formula:

[z^{n%2B1}] H (g(z)) = \frac{1}{(n%2B1)} [w^n] (H' (w) \phi(w)^{n%2B1})

where H can be an arbitrary analytic function, e.g. H(w) = wk.

Lambert W function

The Lambert W function is the function W(z) that is implicitly defined by the equation

W(z) e^{W(z)} = z.\,

We may use the theorem to compute the Taylor series of W(z) at z=0. We take f(w) = w \mathrm{e}^w and a = b = 0. Recognizing that

\frac{\mathrm{d}^n}{\mathrm{d}x^n}\ \mathrm{e}^{\alpha\,x}\,=\,\alpha^n\,\mathrm{e}^{\alpha\,x}

this gives

W(z) = \sum_{n=1}^{\infty} \lim_{w \to 0} \left( \frac{\mathrm{d}^{\,n-1}}{\mathrm{d}w^{\,n-1}}\ \mathrm{e}^{-nw} \right) { \frac{z^n}{n!}}\,=\, \sum_{n=1}^{\infty} (-n)^{n-1}\, \frac{z^n}{n!}=z-z^2%2B\frac{3}{2}z^3-\frac{8}{3}z^4%2BO(z^5).

The radius of convergence of this series is e^{-1} (this example refers to the principal branch of the Lambert function).

A series that converges for larger z (though not for all z) can also be derived by series inversion. The function f(z) = W(e^z) - 1\, satisfies the equation

1 %2B f(z) %2B \ln (1 %2B f(z)) = z.\,

Then z %2B \ln (1 %2B z)\, can be expanded into a power series and inverted. This gives a series for f(z%2B1) = W(e^{z%2B1})-1\,:

W(e^{1%2Bz}) = 1 %2B \frac{z}{2} %2B \frac{z^2}{16} - \frac{z^3}{192} - \frac{z^4}{3072} %2B \frac{13 z^5}{61440} - \frac{47 z^6}{1474560} - \frac{73 z^7}{41287680} %2B \frac{2447 z^8}{1321205760} %2B O(x^9).

W(x)\, can be computed by substituting \ln x - 1\, for z in the above series. For example, substituting -1 for z gives the value of W(1) = 0.567143\,.

Binary trees

Consider the set \mathcal{B} of unlabelled binary trees. An element of \mathcal{B} is either a leaf of size zero, or a root node with two subtrees. Denote by B_n the number of binary trees on n nodes.

Note that removing the root splits a binary tree into two trees of smaller size. This yields the functional equation on the generating function B(z) = \sum_{n=0}^\infty B_n z^n:

B(z) = 1 %2B z B(z)^2.

Now let C(z) = B(z) - 1 and rewrite this equation as follows:

z = \frac{C(z)}{(C(z)%2B1)^2}.

Now apply the theorem with \phi(w) = (w%2B1)^2:

B_n = [z^n] C(z) = \frac{1}{n} [w^{n-1}] (w%2B1)^{2n} = \frac{1}{n} {2n \choose n-1} = \frac{1}{n%2B1} {2n \choose n}.

We conclude that B_n is the Catalan number.

See also

References

  1. ^ M. Abramowitz, I. A. Stegun, ed (1972). "3.6.6. Lagrange's Expansion". Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. New York: Dover. p. 14. http://people.math.sfu.ca/~cbm/aands/page_14.htm. 
  2. ^ Lagrange, Joseph-Louis (1770). "Nouvelle méthode pour résoudre les équations littérales par le moyen des séries". Mémoires de l'Académie Royale des Sciences et Belles-Lettres de Berlin 24: 251–326. http://gdz.sub.uni-goettingen.de/no_cache/dms/load/img/?IDDOC=41070.  (Note: Although Lagrange submitted this article in 1768, it was not published until 1770.)
  3. ^ Bürmann, Hans Heinrich, “Essai de calcul fonctionnaire aux constantes ad-libitum,” submitted in 1796 to the Institut National de France. For a summary of this article, see: Hindenburg, Carl Friedrich, ed (1798). "Versuch einer vereinfachten Analysis; ein Auszug eines Auszuges von Herrn Bürmann [Attempt at a simplified analysis; an extract of an abridgement by Mr. Bürmann]". Archiv der reinen und angewandten Mathematik [Archive of pure and applied mathematics]. 2. Leipzig, Germany: Schäferischen Buchhandlung. pp. 495–499. http://books.google.com/books?id=jj4DAAAAQAAJ&pg=495#v=onepage&q&f=false. 
  4. ^ Bürmann, Hans Heinrich, "Formules du développement, de retour et d'integration," submitted to the Institut National de France. Bürmann's manuscript survives in the archives of the École Nationale des Ponts et Chaussées [National School of Bridges and Roads] in Paris. (See ms. 1715.)
  5. ^ A report on Bürmann's theorem by Joseph-Louis Lagrange and Adrien-Marie Legendre appears in: "Rapport sur deux mémoires d'analyse du professeur Burmann," Mémoires de l'Institut National des Sciences et Arts: Sciences Mathématiques et Physiques, vol. 2, pages 13–17 (1799).

External links